Integrating this model

\[ \frac{dC_{ij}}{dt} = f(C_{ij}, T_{ij}) \]

Where \(C_{ij}\) is sea ice concentration in gridpoint \(ij\), \(T_{ij}\) is the t2m temperature in gridpoint \(ij\) and \(f\) is the polynomial fit on monthly values. \(T\) comes from ERA5.